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The technological performances of metallic compounds are largely influenced by atomic ordering. 
Although there is a general consensus that successful theories of metallic systems should account for 
the quantum nature of the electronic glue, existing non-perturbative high-temperature treatments 
are based on effective classical atomic Hamiltonians. We propose a solution for the above paradox 
and offer a fully quantum mechanical, though approximate, theory that on equal footing deals 
with both electrons and ions. By taking advantage of a coarse grained formulation of the density 
functional theory [Bruno et al., Phys. Rev. B 77, f 55f 08 (2008)] we develop a MonteCarlo technique, 
based on an ab initio Hamiltonian, that allows for the efficient evaluation of finite temperature 
statistical averages. Calculations of the relevant thermodynamic quantities and of the electronic 
structures for CuZn and NisV support that our theory provides an appropriate description of order- 
disorder phase transitions. 

PACS numbers: 



Nowadays ground state properties of metallic alloys 
are routinely calculated by using Density Functional 
(DF) theory which appropriately incorporates elec- 
tronic correlations. However, most finite temperature 
status of art calculations Q are based on classical Ising 
models which allow the accurate evaluation of phase equi- 
libria, but are not able to predict the electronic proper- 
ties. A theory able to cope with both tasks must incor- 
porate quantum mechanics and should be able to explore 
the very large space of alloy configurations [4| in order 
to give reliable statistical averages. Unfortunately, quan- 
tum simulations as the Car-Parrinello Molecular Dynam- 
ics (CPMD) are very hard for the problem at hand and 
the application of such methods the phase equilibria of 
metallic alloys seems beyond the capabilities of existing 
hardware and software. In this Letter we shall outline a 
new approach to the problem, which, similarly to CPMD 
is based on the DF theory and the Born-Oppenheimer ap- 
proximation, but where the relevant quantities shall be 
obtained as statistical averages (by sampling the space 
of alloy configurations) rather than as Boltzmann time 
averages. 

Although the generalization to n-ary systems is 
straightforward, in this Letter we focus on binary metal- 
lic alloys AcBi_c. They shall be studied in the (T, c) 
statistical ensemble defined by the temperature, T, and 
the mean atomic concentration, c. In order to have a 
tractable problem, we shall limit ourselves to the solid 
state and the normal metal regime or, equivalently, to the 
temperatures between the superconducting and the melt- 
ing transitions, Tgc < T < Tm- Furthermore, marten- 
sitic and magnetic phenomena shall not be considered. 

Our first crucial step is to obtain a coarse grained ver- 
sion of the Hohenberg-Kohn DF. If the ions constituting 
the system are considered frozen on the sites of a simple 
lattice at their equilibrium positions, Ri, then the elec- 
trostatic contribution to the DF can be written as a sum 



of local terms plus some bilinear terms involving the local 
moments of the electronic density 0| : 
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In Eq. ([T]), each lattice site is associated with a Voronoi 
polyhedron (VP), i.e. to the set of points closer to 
that particular rather than to any other lattice site, the 
Ui{[pi{r)]) are known functionals of the local electronic 
density, Pi{r), —eqi are the charge multipole moments in 
each VP and the Madelung matrix, M, is determined by 
the crystal geometry only. As commented in Ref. @, the 
coupling between different VP's in Eq. ([T]) is marginal 0], 
in the sense that it has a simple, analytically tractable 
form. Most existing approximations for the exchange- 
correlation energy [8| also consist of sums of local elec- 
tronic density functionals. However, the local kinetic 
contributions to the DF are non trivially entangled by 
the boundary conditions (be) that the wave- function, or 
the Green function Q, must match at the VP surfaces. 
A marginally coupled DF then requires an approxima- 
tion for the kinetic part. The literature reports several 
ways in which such decoupling has been obtained: Ref.fiol 
uses vacuum bc's, while random bc's are employed in 
Ref. lul In our approach the 'true' bc's are replaced 
by appropriate mean field bc's The same procedure 
defines the class of the Generalized Coherent Potential 
Approximations (GCPA), whose prototype, the Isomor- 
phous CPA (ICPA), has been introduce d by Soven 
The ICPA provides an excellent picture [l3| of the spec- 
tral properties of metallic alloys but may lead to incor- 
rect predictions about the electrostatics [l4i. Modern 
GCPA schemes like the Polymorphous CPAjll (PCPA) 
heal the shortcomings of the old ICPA model and pro- 
vide fairly good total energies both for ordered and dis- 
ordered alloy configurations We shall then use the 
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GCPA DF [^: 



(2) 



In Eq. the dependence on the alloy configuration is 
identified by the set of occupation numbers {£,}, with = 
1 (or 0) for sites occupied by an A (or B) atom. In the 
case of a simple lattice, the local functionals w"([pi(r)]) 
parametrically depend on the site chemical occupation 
only, a = A or B. Minimizing Eq. ^ with respect to the 
Pi(r) provides a set of Euler-Lagrange equations coupled 
only through the Madelung potentials, Vi = J^j^ijlj^ 
and the chemical potential fi: 



Suj" 



(3) 



Within the alloy sample specified by {^}, due to 
Eq. ® and to the variational properties of the GCPA 
DF [a^lB, 17 1, the local charge densities for each chem- 
ical component, p"(r), are unique functions [Q] of the 
Madelung potentials. By virtue of the Hohenberg-Kohn 
theorem, the same holds for any other local physical ob- 
servable, 



(4) 



Evidently, also the charge moments, qi, arc unique func- 
tions of the Vi. If, as it can be argued on a physical 
ground, the qV laws, qi — q°'{Vi), are strictly mono- 
tonic &\, then they can be inverted, thus allowing to re- 
cast the first term in Eq. ((21) as a function of the charge 
moments rather than a functional of the full charge den- 
sity, ti;"([pi(r)]) u}"{qi). This gives the desired coarse 
graining. 

Actually, in the metallic state, the qV laws are not 
only monotonic but numerically almost indistinguishable 
from straight lines [1, [3, [l^. Thus, an excellent approx- 
imation for Eq. (m can be obtained by a series expan- 
sion about the zero field (or ICPA) values of the charge 
moments, This gives the Charge Excess Functional 
(CEF) already obtained on a phenomenological ground 
in Ref. [11, 
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CEF 
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where Ui = {(P Cj"^ {q) / dq^) q=q„ takes the values a a or 
as, depending on the site occupation. For a given con- 
figuration, {^}, the minimum value of the functional in 
Eq. (12) or ([5]) corresponds to the total system energy and, 
therefore, provides an ah initio effective atomic Hamilto- 
nian |19j. Unlike most Ising models, such CEF Hamil- 
tonian includes effective interactions at all distances and 
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FIG. 1: Distributions of the Madelung potentials (DMP). 
Calculated DMP's for CuZn (left) and NisV (right) alloys 
at the indicated temperatures. The structures visible in the 
DMP's are associated with the local atomic environments. 
Above the transition temperature (To ~ 4657^ for CuZn and 
To = 1701 ± IIK for NiaV) the DMP's present an overall 
Gaussian shape, while at lower temperatures they resemble 
the T — > limit, consisting of 5-like peaks, one for each non- 
equivalent lattice site. 



n-body interactions up to any value for n [20|. For a 
binary alloy, it is determined by only three parameters 
that can be easily obtained by GCPA-DF calculations. 

In this Letter we implement a coarse grained all elec- 
trons theory based on the GCPA DF and a finite tem- 
perature Monte Carlo (MC) sampling 21, 23] of the al- 
loy configurations space. This theory, in the following 



referred to as GCPA-CEF-MC, consists of three ma- 
jor stages [H]: (i) the GCPA DF is determined by 
T = electronic structure calculations (ii) fo'^ ^ 

given thermodynamic point (T, c) the relevant ensemble 
of configurations is sampled by a MC based on the CEF 
Hamiltonian; for each configuration, the local charges 
and Madelung potentials are obtained by minimizing 
Eq. (15]), whose solutions are numerically indistinguish- 
able [6] from those of Eq. ([2]); (iii) the obtained distribu- 
tion of the Madelung potentials (DMP) (see Fig. [T]) then 
allows for the evaluation of the appropriate ensemble av- 
erages for the electronic properties through Eq. ([4]). As 
an example, we report in Fig. [2] the electronic density 
of states (DOS) for CuZn alloys at several temperatures. 
The stage (n) of our GCPA-CEF-MC provides the en- 
semble averages corresponding to the relevant thermo- 
dynamic quantities and to the properties related with 
the atomic degrees of freedom only. Such is, for instance, 
the atomic Short Range Order (SRO) (see Fig. [31 where 
results for NisV alloys are presented). 

In order to test our GCPA-CEF-MC method, we have 
selected two well studied systems: CuZn 2^ 26, 27 1, al- 
ready discussed by Landau [28*1 as the prototype for his 
theory of second-order phase transitions, and NiaV [2^, 
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FIG. 2: Site resolved electronic DOS's for CuZn alloys. Up- 
per frame: mean DOS contributions from Cu and Zn sites 
at the indicated temperatures. Lower frame: the comparison 
between the mean DOS's from Cu and Zn sites at T = 350K 
(thick black lines) with those from the most extreme chemi- 
cal environments in the GCPA ensemble (ligh lines, the cor- 
responding values for e^Vi are indicated by the labels) shows 
the sensitivity of the GCPA theory to the site environments. 
The lower frame reports also the ICPA site resolved DOS 
evaluated at the same lattice constant (thick gray lines). 



|30|, for which there is competition between different or- 
dered structures. For both alloys thermal expansion has 
been included 3lj. This inclusion, however, implies only 
small quantitative changes to our results. Below Tm, 



both systems present disordered solid solution phases, 
based on the relatively open bcc lattice for CuZn and on 
the close-packed fee for NisV. On decreasing T, both sys- 
tems undergo a transition to the ordered phases, B2 and 
DO22 respectively. The phase transitions have been mon- 
itored by plotting the total energies and the specific heats 
as a function of T (Fig. |4]). The calculations for CuZn 
alloys show the neat evidence of an order-disorder transi- 
tion occurring at Tq ~ with a smooth dependence 
on N, the number of atoms contained in the simulation 
box. In the vicinity of the transition, quite large N values 
are necessary to obtain well converged thermodynamic 
properties for the NisV system. For the largest simula- 
tion box, ordering occurs on cooling between 1692.5 and 
1690 K, while, on heating, the low-T ordered state disor- 
ders between 1710 and 1712.5 K. Our best estimate for 
To in NisV is then 1701 ± 11 K. The transition reveals a 
first-order character, as confirmed by the abrupt change 
in the SRO's displayed in Fig. [3l The sharp first-order 
features of the transformation are gradually suppressed 
by finite size effects in the smaller simulation boxes in- 





FIG. 3: Atomic Short Range Order (SRO) for NisV alloys 
in the ky — plane, at the indicated temperatures. Above 
To, the kx and kz directions are equivalent, due to the cubic 
symmetry of the disordered fee phase. Below To, the cubic 
symmetry is broken, giving rise to the tetragonal symmetry of 
the DO22 structure, with the characteristic peaks at (1,0,0) 
and (1/2, 0, 1). It is remarkable the abrupt change of the 
structure [23I, that on cooling occurs between 1692.5 K and 
1690 K. 



vestigated. 

Our results demonstrate that the GCPA-CEF-MC the- 
ory is able to provide a sensible description of ordering 
phenomena in metallic alloys. Remarkably, we find com- 
plete agreement with the experiment about the low-T 
ordered phases and the theory correctly discriminates be- 
tween first- and second-order transitions. The agreement 
with the experiment about the ordering temperatures, 
however, is not completely satisfactory. The experimen- 
tal ordering temperatures [s^l for CuZn and NiaV, 740 
K and 1313 K, respectively, are in one case lower and 
in the other higher than the calculated values. This is 
probably due to the single-site nature of the GCPA and 
to the approximations 2^ made for the site potentials. 

The methodologies presented in this Letter offer a 
route to alloy thermodynamics that fully includes elec- 
trons and appreciably enlarges the scales at which quan- 
tum mechanics can be applied. The calculations here 
discussed consider up to 2000 atoms and up to 10^^ local 
chemical environments for each point in the (T, c) space. 
It is clear that such impressive figures can be obtained 
only because we have a fixed crystal lattice. The applica- 
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FIG. 4: Electronic contributions to the total energies, E, and 
the specific heats, cv, for CuZn (upper frame) and NiaV 
(lower frame) alloys as a function of the temperature. The 
specific heats (insets) identify the ordering temperatures. To, 
with the respective uncertainties indicated by the shaded re- 
gions, and exhibit the typical A-shape. Different symbols refer 
to different sizes (A'^) of the simulation boxes. CuZn energies 
do not show hysteresis, as expected for the second order tran- 
sition bcc-B2. For the largest simulation box, NisV energies 
show a hysteresis loop, revealing the first-order character of 
the symmetry breaking fcc-D022 transition [2^ . 



tion of our method is then complementary to the CPMD, 
which is very hard to apply to diffusive transitions in the 
solid state. We are confident that, in the next future, 
coarse grained DF theories shall be extended to mag- 
netic metals and to the mesoscopic domain, in order to 
study defects and plasticity. 

The calculations presented have been executed us- 
ing the computational facilities located at the CINECA 
(Bologna, Italy) and at the Consorzio COMETA, project 
PI2S2 (http://www.pi2s2.it). 
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